Primary role of the barely occupied states in the charge density wave formation of 

NbSe 2 
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NbSe2 is a prototypical charge- density- wave (CDW) material, whose mechanism remains myste- 
rious so far. With angle resolved photoemission spectroscopy, we recovered the long-lost nesting 
condition over a large broken-honeycomb region in the Brillouin zone, which consists of six saddle 
band point regions with high density of states (DOS), and large regions away from Fermi surfaces 
with negligible DOS at the Fermi energy. We show that the major contributions to the CDW for- 
mation come from these barely occupied states rather than the saddle band points. Our findings 
not only resolve a long standing puzzle, but also overthrow the conventional wisdom that CDW is 
dominated by regions with high DOS. 
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2H structured niobium diselenide (NbSe2) is one of 
the most studied materials for its prototypical supercon- 
ductivity (SC) and two dimensional charge density wave 
(CDW) 0, S, & i, i, i]. Based on experiments con- 
ducted on NbSe2, understanding on some of the most 
basic properties of SC and CDW was reached, such 
as the internal electronic structures of the vortex cores 
7 1 .the anisotropic s-wave 0, @, 0] and multi-band SC 
S , i lfj, llll , the anisotropic electron-phonon interaction 
6|, |l2| , and the phonon softening at the CDW transition 



13j . Nevertheless, the very mechanism of the CDW in 
NbSe2 itself has been mysterious and controversial for 
over three decades 0, [l2 , 14, Tj^, [T|| , even though it is 
one of the very first two dimensional CDW materials dis- 
covered (l7j . 

In a typical picture of CDW, there are usually paral- 
lel Fermi surface (FS) sections, which guarantee a large 
number of electrons could be scattered from one parallel 
side to the other by softened phonons with a fixed mo- 
mentum transfer corresponding to the ordering wavevec- 
tor(s). As a result, the charge susceptibility at the order- 
ing wavevector(s) will be enhanced, and CDW instabil- 
ity could be induced in collaboration with lattice. This 
gives the so called nesting condition of CDW. Similarly, 
saddle band points, where the singularities of density of 
states (DOS) are located, were proposed to cause CDW 
as well Recently, the authors discussed another so- 
called Fermi-patch CDW mechanism for some polaronic 
systems, where large DOS at Ep hovers over an extended 
area of the Brillouin zone (BZ ) and work collectively to 
favor certain CDW order [l8(. However, the CDW in 
NbSe2 cannot be explained by any of these scenarios: nei- 
ther can its ordering wavevectors match the FS sections 
0, EE E3|) nor match the saddle band points 15, Hij |. 
even though some experiments showed that the CDW 
gap does open near the saddle band points [T^] ; and there 



are no Fermi patches. The latest experiment identified 
the CDW gap in the hashed circles of Fig.[4ja), some of 
which fulfill the nesting condition 12j. However, there 
are still some regions that do not. The inapplicability of 
these common models for CDW in such a classic com- 
pound suggests that our current understanding of CDW 
need to be amended. 

All the existing CDW mechanisms are centered on the 
regions with high DOS at Ep in the BZ. In this Letter, 
we show that the CDW related spectral weight suppres- 
sion with large energy scale occurs in a large broken- 
honeycomb region mainly away from the FS's, including 
a large momentum area where the DOS is very small at 
Ep. Recovery of the region makes the nesting conditions 
naturally fulfilled, which explains the previous controver- 
sies. Particularly, we prove that it is the states in these 
low-DOS regions that play a primary role in the CDW 
formation. Such an unexpected finding overthrows the 
conventional wisdom that CDW is dominated by the re- 
gions with high DOS. The electronic states in the entire 
BZ have to be examined no matter how low the DOS 
they possess. 

High quality NbSe2 single crystals were synthesized by 
the vapor-transport technique. The data were collected 
with 21.2eV photons from Helium discharge lamp and 
BL9 of Hiroshima synchrotron radiation center by the 
Scienta R4000 electron analyzers. The angular resolution 
is 0.3° and the total energy resolution is ~ lOmeV. In 
the temperature dependence studies, measurements were 
conducted in a cyclic way to guarantee no aging effects in 
the spectra. All experiments were performed in the ultra- 
high vacuum (better than 3xl0 -11 mbar in the helium 
lamp system and ~ 5 x 10~ n mbar at BL-9) and in a 
short time. 

Photoemission intensity maps integrated within 
±10meV around the Fermi energy (Ep) at the normal 
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FIG. 1: (Color online) (a) Photoemission intensity maps (in- 
tegrated within 20meV around Ef) at the CDW state (15K) 
and normal state (45K) are compared. The image taken at 
15K is 6-fold symmetrized, (b) The photoemission intensity 
along Cut 1 marked in panel (a), (c)-(e) The typical ARPES 
spectra for NbSe2 taken at normal state (45K) along the cor- 
responding cuts. The thicker curves are the spectra at the 
numbered momenta in panel (a), and the open circles repre- 
sent the location of the flat band. 



(45K) and CDW states (15K) are compared in Fig.QJa). 
Three FS pockets, one hexagon around T and two around 
K point can be directly distinguished. Furthermore, be- 
cause of the coupling between two NbSe2 layers in a 
unit cell, two split bands can be identified [see Fig.QJb)], 
which gives double- walled FS's [3, H3, E|. Although 
there are long straight sections in the T pockets, none 
of them could be connected by the CDW wavevector 



Q,-~ai*/3=0.688 A (1=1,2,3) [22]. 

Figs.[ljc-e) illustrate the normal state spectra of 
NbSe2 along three cuts across FS's. None of them ex- 
hibits a sharp quasiparticle-like lineshape, even at the 
Fermi crossings. However, they do show clear disper- 
sions, which make the FS's well defined. These have been 
found to be typical signatures of polaronic systems, such 
as Lai.2Sri.sMn207 and K0.3M0O3, where the weight of 
the quasiparticle peak is vanishingly small, and its disper- 
sion is renormalized to the vicinity of the FS. However, 
the majority of states are distributed to much higher 



energies as incoherent spectral weight by scattering of 
multi-phonons[2^, 24 1, which interestingly, still follows 
the bare band dispersion. We note that in NbSe 2 the 
spectral weight at Ep quickly drops when away from 
the FS's, unlike a sibling compound 2H-Na2;TaS2, where 
large spectral weight at Ep is observed almost over the 
entire BZ [l8j]. Moreover, midway between T and K, 
where the saddle band points stand, we indeed observe 
the flat dispersion as marked by the open circles in 
FigS-QJc) and [He). Even though the spectral centroid 
is below Ep, the DOS is quite high here. 

The intensity maps do not show any typical CDW ef- 
fects, such as the folding of the FS's. However, through 
careful examination of the spectra, the suppression upon 
CDW is observed in various regions of the BZ. Several 
representative spectra are plotted in Figs.[2ja-d) for il- 
lustration. The corresponding momenta are marked in 
Fig.QJa) and their corresponding spectra have been high- 
lighted in Figs.QJd-e). The left panel of Fig. [2a) shows 
a spectral lineshape evolution as a function of tempera- 
ture near the saddle band point. The spectra taken at 
45K and 35K, which are both above the CDW transition 
temperature 33K, nearly coincide with each. Once enter- 
ing the CDW state, the spectrum is suppressed notice- 
ably over about 80 meV below Ep while the leading edge 
shifts negligibly. This suppression energy scale is much 
higher than the ksT c of 3 meV, which is the energy scale 
of CDW effects on the quasiparticles. Such low energy 
effects might be passed to the incoherent spectral weight 
at high energies by multi-phonon scattering processes as 
observed here. 

The suppression at saddle point agrees with previ- 
ous observations [13]. However, the main discovery here 
is that such CDW related suppression also occurs over 
large momentum regions inside the T pocket [Fig.f2Jb)] 
and around the M point [Fig. HJc)], though the spectral 
weights at Ep are low there [see Figs.QId-e)]. Whereas, 
for the spectra taken inside the K pocket [Fig. EJd)], no 
such suppression is observed. To show the suppression is 
not due to thermal broadening, all spectra are divided by 
the resolution convoluted Fermi functions at correspond- 
ing temperatures and then multiplied by that of 10K, 
before they are compared in the right panels of Figs. [2]^ a- 
d). While the suppression at momenta No. 1-3 is ob- 
served upon CDW formation, it is absent at momentum 
No. 4. This is further illustrated by the spectral weight 
suppression ratios in Fig.[2]Je). We stress that such sub- 
tle spectral weight suppression is assured to be purely 
dependent of temperature instead of the aging effects or 
noises. Besides the high sampling points and fine statis- 
tic, at the end of the temperature evolution study, the 
first spectrum is always repeated to make sure they over- 
lap, as testified by the two coincident 15K spectra in 
Fig.[2a). 

Following this analysis, more than 1,200 pairs of the 
normal and CDW state spectra were compared one by 
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FIG. 2: (Color online) (a)-(d) Typical temperature dependence of the original (left) and temperature broadening effects removed 
(right) spectra in the normal and CDW states at momenta No. 1-4 marked in Fig.[TJa) respectively (see text for details), (e) The 
spectral weight suppression ratios (compared to spectra taken at 45K and over (EF-80meV , Ef)) for different temperatures 
in fig (a-d) respectively. 



one over more than 1/6 of the BZ. Figs.^c-h) give some 
more examples, where the left panels show the raw spec- 
tra, and right panels show the corresponding spectra with 
thermal broadening removed as practised in Fig. 2. Even- 
tually, the momenta regions where the spectral weight is 
suppressed can be determined accurately, as marked in 
Fig-EJa) by bars. These regions are mainly composed 
of the large domains inside the T pockets, around the 
M points, and around the saddle band points. Such 
distribution has been confirmed by many samples, and 
the results based on two of them are shown here. Note 
this distribution includes some "barely occupied" regions 
where the DOS at Ep is very weak. We emphasize that 
the weak spectral weight and its suppression are not due 
to the umklapp scattering of the CDW order, since it is 
not observed in the FS or dispersion. 

To further visualize the distribution of the weight sup- 
pression, the intensity map at 15K [integrated over (Ep- 
80meV, Ep)] is subtracted from that at 45K, and the 
result is shown in the upper half of Fig.[3^b). Note the 
thermal broadening effects have been removed from the 
spectra following the standard symmetrization method 
before the processing [25| . After normalized by the inten- 
sity map at 45K, the resulting spectral weight suppres- 
sion ratio map [lower part of Fig.[3^b)] further highlights 
the weight-suppressed regions, which are consistent with 
those determined in Fig.[3]Ja). Note we have surveyed 
the entire BZ, and only observed suppression instead of 
spectral weight gain in the measured energy range. Con- 
sidering the conservation of the total electrons, the sup- 
pressed weight should be shifted to some higher binding 
energies. In this way, the electronic energy gain for the 
CDW transition is naturally fulfilled. 

The distribution of the CDW induced spectral weight 
suppression makes a broken-honeycomb-shaped structure 



in the BZ, as shown by Fig.[4ja). The suppression over 
the saddle band regions (surrounded by the dotted lines) 
only account for 43% of the total one. In other words, al- 
though the DOS is higher around the saddle band points, 
it only plays a secondary role in the CDW formation 
due to its relatively small momentum space. Moreover, 
the broken-honeycomb regions can be well connected by 
the CDW wavevector Qi's [represented by the double- 
headed arrows in Fig. 4(a)]. In previous studies, since 
only a small portion of the weight-suppressed region was 
discovered (hashed circles) , multiple scattering has to be 
invoked to barely make some of the regions connected by 
the two dashed arrows in Fig.HJa) [12| . However, this 
multiple scattering assumption is not needed, once the 
large weight-suppressed region is recovered here. 

To further demonstrate the contribution to the CDW 
instability of the broken-honeycomb-shaped region, we 
calculated the weighted joint DOS, 

C(q, uj)= M (k, k + q)A(k, u>)A(k + q, u)dk 



where M = { 



1 (k + q, k € gapped region) 



It has 



(otherwise) 

been shown before that C (q, 0) could effectively describe 
the phase space for scattering of electrons from the state 
at k to the state at k+q by certain modes with wavevector 
q, and it is closely associated with the electronic suscep- 
tibility 0I2& HI- The matrix element M(k, k + q) is in- 
troduced assuming that only the states that directly cou- 
ple (thus suppressed) with the softening phonon branch 
around Qi [13j , would contribute to the CDW instability. 
C(q,0)'s computed on the honeycomb regions [Fig.|4jb)] 
and on the saddle band regions [Fig.[4jc)] are plotted in 
the same scale. Clearly, the peaks at Qi's are very promi- 
nent in Fig.|3Jb), whereas they are not so pronounced in 
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FIG. 3: (Color online) (a) The illustration of the regions 
where gap opening occurs upon CDW transition in the first 
BZ. The gapped regions determined by comparison of the 
spectra below and above Tcdw one by one. (b) The spectral 
weight suppression map and its corresponding ratio map (see 
text for detail) are false color plotted, where the dashed lines 
contour the obviously suppressed regions and the black ones 
indicate the FS's. (c)-(h) Typical spectra (left) and the cor- 
responding thermal broadening removed ones (right) in the 
normal and CDW states at momenta marked in panel (a). 




FIG. 4: (Color online) (a) The distribution of CDW related 
weight suppression over the entire BZ sketched from the spec- 
tral weight suppression map, which are marked by the or- 
ange broken-honeycomb-shaped structure. The gapped re- 
gions could be well connected by the CDW wavevectors, Qi's 
(i— 1,2,3), as indicated by the double-headed arrows. The 
gapped regions proposed by Kiss et al. and around the sad- 
dle band points are marked by the hashed circles and the 
dotted line surrounded rectangles, (b)-(c) Two dimensional 
joint DOS results on the honeycomb regions and the saddle 
band point regions, respectively, (d) Comparison of the cor- 
responding joint DOS results along the dashed lines in panels 
(b) and (c). 



Fig-BJc). This is more clearly illustrated in Fig.[4jd), 
where the C(g,0)'s along the T-M direction are com- 
pared. The saddle band regions only contribute about 
17% of the total C(Qi,0) computed over the broken- 
honeycomb region respectively [281 ]. Therefore, the nest- 
ing of states with low DOS at Ep indeed dominates the 
charge instability at Qj's here. 

The data here also provide some clues on the driving 
force behind the uncoventional CDW transition in NbSe2 . 
For the conventional Peierls transition at one dimension, 
the CDW is driven by the electronic energy gain, and it 
is reflected by a peak in the charge susceptibility [29| . 
Although it might be more complicated in NbSe2, the 
observed peak in the weighted joint DOS does hint that 
the electronic energy gain might be the dominating factor 
here as well. 

In conclusion, we have found that the CDW re- 
lated spectral weight suppression in the large broken- 
honeycomb region mainly away from the FS's for NbSe2, 
which naturally fulfills the nesting condition, and resolves 
a long-standing puzzle of this classic CDW material. De- 
spite of the weak DOS in most of this region, their pri- 
mary role in the CDW formation is demonstrated. Our 
findings show that collective contribution of the weakly 



occupied states can play the primary role in CDW for- 
mation in some cases, and thus should not be overlooked, 
which amends the conventional wisdom of CDW. 
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